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We consider the J/^f production in proton (deuteron) - nucleus collisions at high 
energies. We argue that the production mechanism in this case is different from 
that in pp collisions due to gluon saturation in the nucleus and formation of the 
Color Glass Condensate. At forward rapidities (in the proton fragmentation region), 
the production of J/^f is increasingly suppressed both as a function of rapidity and 
centrality. On the other hand, at backward rapidities at RHIC (in the fragmenta- 
tion region of the nucleus) the coherent effects lead to a modest enhancement of the 
production cross section, with the nuclear modification factor Rj/q, increasing with 
centrality. We find that the J/^> production cross section exhibits at forward rapidi- 
ties the limiting fragmentation scaling established previously for soft processes; in 
the energy range studied experimentally, it manifests itself as an approximate "xf 
scaling" . 
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I. INTRODUCTION 

Understanding the J/ty production mechanism is one of the challenges of QCD. On 
one hand the charm quark mass is quite large on the typical QCD scale of Aqcd, which 
makes the use of perturbative QCD meaningful [l] since the long distance dynamics is 
effectively decoupled 3|. However the size of this system and the inverse of the binding 
energy are not small enough to suppress significant non-perturbative contributions. Indeed, 
perturbative QCD fails in describing the differential J/ty production cross section and the 
polarization. Different mechanisms were suggested to explain the existing experimental data. 
Unfortunately all of them so far have encountered problems in describing at least some of 
the observables (for a recent review, see jsj). In the context of high energy nuclear physics, it 
is important to understand well the mechanism of J/ty production also in nuclear processes 
since J/^f suppression in heavy- ion collisions could serve as a signal of the Quark-Gluon 
Plasma formation ^j]. 

One of the long-standing puzzles is the lack of x^ scaling of the nuclear modification factor 
{xi is the Bjorken variable corresponding to the nuclear target parton distribution) in J/ijj 
production off nuclei. Even though this scaling is expected to hold in the parton model, the 
data from CERN H and FNAL 6] fixed target experiments are in violent contradiction with 



this expectation. The absence of %2 scaling has become even more dramatic at RHIC |7j]. 
Instead of the badly broken x 2 scaling, the data instead exhibit an unexpected approximate 
scaling in the Feynman xf variable. 

It was realized long time ago [8] that this lack of %2 scaling, and thus the violation 
of QCD factorization, is caused by multi-parton (higher twist) interactions in the nuclear 
target. Several specific mechanisms of this type were considered over the years [9|, HO, 
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In this paper we would like to re- visit the problem of J production in proton - nucleus 
collisions at high energy basing on the novel Color Glass Condensate picture of the nuclear 
wave function at small x. In this approach, the strength of the color field inside the nucleus 
is proportional to the saturation scale Q 2 s (x2) determined by the density of partons in the 
transverse plane. It is a growing function of the collision energy and the atomic number of the 
nucleus. Experimental data indicate that at RHIC kinematics Q s ^> Aqcd which implies 
that the inter-nucleon interactions play a little role in pA interactions at high energies. 
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Therefore, at high energies a nuclear color fielc 
(process independent) dimensional scale Q s (x 2 ) 



can be described by only one universal 
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quarks in this framework has been previously considered in several papers 
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In the previous publications [2CJ, |21[ , we have argued that there are two different dynam- 
ical regimes of heavy quark production at high energies depending on the relation between 
the saturation scale and the quark's mass m. When Q s Cm the heavy quark production is 
incoherent, meaning that it is produced in a single sub-collision of a proton with a nucleon. 
This case can be treated within a conventional perturbative approach. In the opposite limit 
of Q s ^> m the heavy quark production is coherent since the whole nucleus takes part in 
the process. In this case the heavy quark production is sensitive to a strong color field 
(CGC) which violates the decoupling of the subprocess of heavy quark production from the 



dynamics of partons in the nuclear wave function 20, 26]. The reason is that the decoupling 
theorems can be applied only when the heavy quark mass is much larger than the typical 
hadronic scale, which is of the order of Q s {x) at high energies. 

The goal of this paper is to address the problem of J/^J production at high energies. It 
is organized as follows. In section |H] we argue that at high energies the time of interaction 
of the projectile proton with the target nucleus at rest is much smaller than the time of 
heavy quark pair production and subsequent formation of a bound state. This will allow 
us to use the eikonal approximation and derive in section ITTT1 the cross section for the J/\I/ 
production in pA collisions at high energies. In section HVl we study the derived expression in 
two different kinematical regions. We show that at Q s <C M$ the J/ty the cross section is an 
increasing function of centrality since the dominant contribution comes from the two-gluon 
exchange process. At Q s 3> M$ multiple re-scattering as well as quantum evolution lead to 
the suppression of J/^f production both as a function of energy/rapidity and centrality. We 
also point out that the approximate xp scaling observed in SPS and FNAL data emerges 
naturally in our approach. However, it is seen to be a consequence of the slow dependence 
of the gluon distribution on energy, and so is broken at higher energies. However, as the 
energy increases (e.g. to the LHC energy range) the xp scaling is restored owing to the onset 
of gluon saturation in the incident proton. We compare our calculations with the available 



experimental data in Sec. |V| and conclude in Sec. EI 
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II. DIFFERENT PRODUCTION MECHANISMS OF J/¥ ON NUCLEUS 

A. Relevant time scales 

Consider charmonium production in a pA collision. A cc pair is produced over the time 
/ ( 2777. ) in its center-of-mass frame. In the nucleus rest frame this time is Lorentz-dilated 
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1 Eg 

TP ~ (2m c ) (2m c )' 1 j 

where .E^ is a parent gluon's energy. This time scale should be compared to the typical 
interaction time T mt ~ Ra/c. Eq. (0) shows that at very high energies the production time 
of cc pair can be much larger than the interaction time Tp ^> T mt . This a general property 
of all hard processes at high energies: they develop over a long time Tp * |27| . 

This formula can be rewritten in terms of the Bjorken variable associated with nucleus, 
X2- Note, that the gluon takes fraction x\ of the proton's energy E g = X\E V . Also, by 
four- momentum conservation, (2 m c ) 2 = xix 2 s = 2x\x 2 E p , where is the nucleon 
mass and E p is the proton's energy. Thus, it follows from (0) that 

T P « . (2) 

2M N x 2 V ; 

At RHIC, in the center-of-mass frame x 2 = (m c / \/s) e~ y = 6.5 • 10~ 3 e~ y , where we intro- 
duced rapidity y. Therefore, 

t p (RHIC) w 15e y fm. (3) 

Equation (J3J) implies that at forward rapidities y > 1 one can indeed assume that the 
proton interacts coherently with the whole nucleus (similar estimates for the fixed target 
energies can be found in |2£j). In this case the transverse size of the cc pair is fixed during its 
propagation through the nucleus and we can apply the eikonal approximation for calculation 
of the scattering amplitude 

Finally, the J/\l/ wave function is formed from the initial cc pair. This lasts ~ 2l{M^i — 
My) in the J/\l/ rest frame. In the nucleus rest frame the J/^I formation time is |29j, |30j 

2 E g 



Sometimes one introduces the "coherence length" l c — Tp/c. 
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Since the J/^f bounding energy is much less than its mass, the J/^> production time is much 
shorter than its formation time: rp < r^. A more accurate evaluation of the formation 
time can be performed with the help of the spectral representation for the J/^f propagator 
by using the experimental data on e + e~ annihilation into charm quarks this leads to 
the J/ 1 ® proper formation time of 0.45 fm. Thus, at RHIC 

Tf(RHIC) ss 41e y fm; (5) 

the J/ 1 ® wave function is therefore formed outside of the nucleus at rapidities y > — 2. 

The relationships between the three relevant time scales r int , Tp and Tp depend on the 
collision energy \fs and the rapidity y. In the next subsection we classify all possible 
situations in the RHIC kinematical region (y/s = 200 GeV). 

B. Coherent versus incoherent J/^f production 

1. Forward rapidities 

At (pseudo)-rapidities y > 1 rp ^> Tp ^> r int . This implies that the cc pair scatters 
coherently off all the nucleons along its trajectory in the nucleus, see Fig. [0 Therefore, 
the process of J/^f formation proceeds through the following three stages in the nucleus 
rest frame. First, way before the collision with the nucleus, the fast proton develops a 
cloud of virtual partons which includes one cc pair (in the leading order in a s ). In the 
light-cone perturbation theory this is described by the valence quark and the virtual gluon 
wave functions which are explicitly displayed in the next section, see IIII Bl Second, the 
coherence of the cloud is destroyed by interaction with nucleons in nucleus. Due to the large 
production time Tp the scattering matrix can be diagonalized in the color dipole basis: the 
transverse size of the cc is fixed during the interaction. We calculate the ccA amplitude in 
section UlI Al Third, J/^> is formed far away from the nuclear remnants. No nuclear effects 
are expected at this stage. 

2. Central rapidities 

At rapidities —1 < y < the production time becomes smaller than the nuclear size which 
implies that the cc pair scatters coherently off a few nucleons. This leads to the enhancement 



6 




of J/ 1 ® production since in that case the main contribution to the scattering amplitude arises 
from the diagram shown in Fig. a) which is enhanced by an additional power of A 1 / 3 with 
respect to the diagram Fig. [2Kb) which describes J/*f> production in pp collisions. 

The diagram Fig. Efa), where the J/\I/ is produced by double gluon exchange is paramet- 
rically enhanced compared to the diagram shown in Fig. Efb) where the J/^f is produced 
by one gluon exchange. Indeed, let us for a moment concentrate on J/ty production in a 
quasi-classical approximation where the coupling is small a s (Q 2 s ) <C 1 and together with 
atomic number A ^> 1 it forms a resummation parameter a 2 A 1 / 3 ~ 1 j^. The diagram (b) 
in Fig. El is of the order a\A 1 ^ 3 ~ a 3 while the diagram (a) is of the order of a & s A 2 ^ 3 ~ a 2 .. 
Therefore, the diagram (a) is enhanced provided that the nucleus is sufficiently large. 

This conclusion remains valid beyond the quasi-classical approximation. The gluons 
emitted in the course of quantum evolution get resummed into gluon distribution functions, 
which therefore grow fast as x decreases. The resummation parameter a s xG(x, Q 2 ) A 1 / 3 
becomes large even for the proton [161 ] . At small enough x the diagram (a) dominates the 
J/ty production even for the scattering off light nuclei A ~ 1. The general effect of the 
enhancement of double-gluon exchange diagrams in hard processes on heavy nuclei has been 
already pointed out by in (^J [3] • 

Although the dipole model qualitatively describes the effect of the enhancement of J/ty 
production at — 1 < y < it cannot be applied at y < due to the effect of finite production 
time. In other words to get a reasonable description of the process one has to include the 
absorption corrections of the color dipole cc in the nuclear medium (note that rp is still 
much larger than Ra)- In the present publication we are going to analyze only y > region 
at RHIC. However when comparing to the experimental data from CERN and FNAL we 
will have to correct the results by the absorption factors, see Sec. IV Bl 



(a) (b) 

FIG. 2: Two mechanisms of J/^f production: (a) by two-gluon exchange; (b) by one-gluon exchange 
and one-gluon emission. 



To quantify the effect of the finite production time and to specify the region of applica- 
bility of eikonal approximation, we can consider the longitudinal nuclear form factor F^(q z ), 
which takes into account the quantum interference in the longitudinal direction at a finite 
longitudinal momentum transfer j^l (this formfactor is assumed to equal unity in the dipole 
model). It is defined as 

2 



//*oo 
d 2 b / dzp(b,z)e iq * z 
J — oo 



(6) 



where q z is the longitudinal momentum transfer and p(b, z) is the density of nucleons in a 
nucleus. The typical value of the longitudinal coordinate in (JBJ) is z ~ Ra- Therefore, the 
integral over z vanishes due to rapid oscillations of the exponential factor unless q z < 1/ 'Ra- 
On the other hand q z ~ l/ip. Thus, Fa^z) is sizable (it is normalized to unity) when 
t p > Ra and is small otherwise. 

Assuming for simplicity a Gaussian parameterization of the nuclear density we obtain 
^Klz) — e~ R A q z/ 3 3^. In the Fig. El we plot the longitudinal form factor for different 
energies. We learn from Fig. Elthat the corrections to the dipole model due to a finite Tp are 
of the order of 10% at y = at RHIC (y/s = 200 GeV), but already at y = — 1 they are as 
large as 50%. Therefore, the dipole model should provide an accurate result for y > which 
correspond to xp > 0. At Fermilab fixed target experiments (y/s = 38 GeV) the dipole 
model is applicable (corrections < 10%) at rapidities larger than ~ ln(200/38) = 1.7 which 
correspond to xp ~ 0.4. At SPS (y/s = 19 GeV) the dipole model has to be corrected even 
in the most forward region. 
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FIG. 3: The longitudinal form factor F\ as a function of rapidity y for RHIC, Tevatron and SPS 
energies. 

3. Backward rapidities 

At rapidities y < — 2 the coherence is completely lost and the process of J/\l/ production 
in pA collisions becomes similar to the one in pp collisions. All of the dependence on A arises 
then from the propagation of the produced cc pair and J/\l/ through the nuclear matter. 

III. CROSS SECTION FOR J/* PRODUCTION AT FORWARD RAPIDITIES. 
A. Propagation of the cc pair through the nucleus 

The first step in calculation of the J/\l/ production cross section, see Fig. ^ , is the 
calculation of the scattering amplitude of q v cc state off the nucleus with the projection on 
the color neutral state after the cc passes the nucleus. In our analysis we will assume, for 
the sake of simplicity, that the valence quark is a spectator. 

It is straightforward, although laborious, to calculate all possible gluon attachments to 
quark and antiquark in Fig. ^ In the large N c approximation the cc amplitude for a single 
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scattering (Fig. |2J) projected on the color singlet state reads 



M^x^yJ = [{x-yf + (x -y ) 2 ] . (7) 



The saturation scale Q 2 S in (fTU|) is given by 
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Ql{x) = ^^ P T{b)xG{xA/x\ (8) 

where the gluon distribution function in the nucleon reads 

xG (xA/x 2 ) = ^^\n^, (9) 
7i xrn z 

where \x is some infrared cutoff. For a spherical sharp-edge nucleus T(b) = 2\J R 2 — b 2 . 

Assuming that the scatterings off the individual nucleons are independent, one can gen- 
eralize (0) to include multiple rescatterings (see Fig. |TJ) in a straightforward way. The qqA 
amplitude then takes the form 

M&ZoiV'Vo) = 1 " ex p{ ~\{t) 2 [^-yf + feo-^o) 2 ] 2 } ( 10 ) 
We checked that the color factor exponentiates using the FeynCalc package j^sl ]. 

B. Wave functions 

In view of the arguments given in the previous section, at forward rapidities the processes 
of formation of the proton wave function, the subsequent rescatterings of the proton's partons 
inside the target nucleus and the formation of the bound cc state are well separated in time 
in the nucleus rest frame. Hence, to proceed we need to know the light-cone wave functions 
of the valence quark, the virtual gluon and the J/^f. In the A+. = light cone gauge the 
valence quark's wave function in the configuration space reads 23] 

= / ^e- iq -%^Ug) =9T a ^^, (11) 



Averaging the square of Eq. (jlljl over the quantum numbers of the initial quark and summing 
over the quantum numbers of the final quark and gluon we obtain the familiar gluon radiation 
kernel of a dipole model 



$qv-+qv9\£.li22) oy / . ' —'/'■■'/ ' —>. ' ~ <i ,■—</, — -' ' .2-2 
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where z_ x and z_ 2 are the transverse coordinates of the gluon in the amplitude and in the 
complex conjugated amplitude correspondingly, see Fig. [TJ 



The light-cone wave function of a virtual gluon of momentum q reads j21j 



Vg*->qq(k,k-q,a£) = T 75- =■ (S r y (fc - a q) ■ e [r(l - 2a) + A] +r6 r /m(l + rA)), 

(13) 

where A; is the produced quark's transverse momentum, m its mass, a = k + /q + is the 
fraction of the gluon's light-cone momentum q + it carries, r and r' are the quark and the 
antiquark helicities correspondingly. Projecting it onto the vector meson wave function 
and summing over the polarization and helicity indices one can find the overlap function 
*&v(kj a ) * ^g*{k, a ) i n momentum space j^lj which can be Fourier transformed to config- 
uration space. In the non-relativistic approximation \Ev(r, a) oc 5(r) S(a — \) the overlap 
function in the configuration space takes the form [42] 



1, l3T ee M^ m 3 r 2 



where r = x — x , is the modified Bessel function, T ee ~ 5.26 KeV is the leptonic width 
of J/ty and we will include two delta functions directly into the expression for the cross 
section. 

C. Cross section 

Using (|12j) . (|14|) . and (jlOj) we can obtain the inclusive J/^> production cross section 

x ^ * (x-Xq, a) % * V*. (y - y Q , a) M (x, x y, e -fc-(a-2) 
x6((xq- z x ) + a(x - x )) 5((^ - z 2 ) + a{y - yj) 

x<5((x -u) + oi(x-x J0 ))8((y J0 - v) + a(y-y J0 )) (15) 

where u and y_ are the J/^ coordinates in the amplitude and in the complex-conjugated one 
respectively. Performing the integration over u and v using the last two delta-functions in 
(fT5Jl somewhat simplifies this expression: 

da f j2i f j2 r2 ol s C f z x -z 2 fd 2 xd 2 y f d 2 x d 2 y f 

MTy = J d " J d Zld Z2 ~ SS J W J ~WTJ da 
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x^*^ g .(x-x ,a)%*^* g .(y-y ,a) M{x,x y_,y o )e 

xS((x - ?i) + a(x- x )) 5((y o - z 2 ) + a(y- y Q )) (16) 

Equation (jlfij) gives the differential cross section for J/^> production in pA collisions in 
a quasi-classical approximation to the nuclear color field, at large N c limit and neglecting 
relativistic effects in the J/\I/ wave function. The inclusion of high energy quantum evolution 
effects is important for phenomenological applications at RHIC, but is a quite difficult 
problem. Fortunately, the total inelastic cross section is not very sensitive to the evolution 
effects as we argue below. Another important reason to focus on the total cross section is 
that it is much less dependent on a model which we choose to describe the vector meson 
wave function. 

The total J/^f production cross section per unit rapidity is found by integration over the 
transverse momentum p in ()16j) . It yields the delta function (2ii) 2 5(z l — z 2 ). It is convenient 
to introduce the cc-dipole transverse separation two- vector in the amplitude r = x — x and 
the complex-conjugate one r' = y — y Q . Now, upon substitution of M. into (|TB|) we get 

d a f j2i f j2 a s Cf 1 f d 2 r d 2 r' f d 2 Xod 2 y 



T y = J db J dz —?i^?i^r l(U ' 



x ^ * # 9 *(r, a) % * ^.(r', a) (l - e^(%) [feo-^) 2 +(r-r'+&-i ) 2 ] 2 
x <K^o — z + ar) 5(y Q — z + ar') (17) 

1 _ e -3(%)V 2 +(i-") 2 )(r-r') 2 ] 2N ) ( 18 ) 



Using an explicit formula for the overlap function (fT4*|) in ()18|) we evaluate r[ integral in the 
approximation r' <C r. The result is 

drr 3F f°° / (Q*i*2K\\ 

= S A xG(x 1 ,Q 2 ) 2 ee / dCC 5 K 2 (C)(l-e^^r } ). (19) 

where we used © and introduced a dimensionless variable ( = mr. The gluon distribution 
in the proton xG is evaluated at the scale Q 2 = M 2 + kA 1 ^ 3 , with k a parameter to be 
fixed by experimental data. For simplicity we assumed that the nucleus profile function is 
T(b) = 9(Ra — b) where Ra is an effective, centrality-dependent radius determined by the 
Glauber analysis of pA and dA interactions, see e.g. 44]. 
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IV. INTERPLAY OF TWO SCALES: Afy AND Q s . 

A. The effect of quantum evolution 

Equation (fTHjl gives the desired result for the total cross section of J/^/ production in 
the quasi-classical approximation. In this approximation the saturation scale Q s is given 
by (JHJ). As the result of quantum evolution the saturation scale acquires its energy /rapidity 
dependence: 

Q^) = A^i/3 (^)%\ (20) 

where y is the rapidity measured in the center-of-mass frame. The value of Q 2 S thus increases 
from the initial value A 2 A 1 / 3 given by (jHJ)- The value of A = 0.63 GeV is fixed by DIS 
data. The rate of increase is set by the factor A. It is constant in the leading logarithmic 
approximation A = Aa s N c /-K [45]. Various inclusive quantities at small x at RHIC and 



HERA are well fitted with A « 0.25 



obtains from the Renormalization Group improved 5 1| BFKL equation ,52] . Moreover, the 
value of A is approximately constant in the relevant for us range of virtualities, so we assume 
that A(M^) « 0.25. 

In our discussion we assume that evolution of the gluon density in the proton is linear 
since the saturation scale in the proton Q 2 is A 1 / 3 times smaller than in the nucleus. This 
is a justified approximation in the RHIC kinematical region where the values of X\ are such 
that the proton wave function is dominated by the Fock states with a relatively small number 
of gluons. However, at the LHC the gluons in the proton will also likely be saturated. We 
will discuss the implications of this in Sec. II V CI 



This is close to the value one 



B. Enhancement versus suppression of J/^f production 

Let us now investigate the J/^> production cross section in two kinematical regions: (i) 
M^j > Q s and (ii) < Q s . The nuclear effects are usually expressed in terms of the nuclear 
modification factor defined as 
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The J/^f production cross section in pp collisions can be obtained by expanding the exponent 
in (HH 



dy p y (27r) 2 a em M^ V M$ J 

where Q sp is given by (f2T)|) with A — 1. 

In the region (i) we can expand the exponent in (fTTH) and using ()22j) we find 

Rpa(J/V) = A 1 / 3 ~ , > Q s , (23) 

which means that the J/^f production is enhanced at backward rapidities at RHIC and it 
is stronger for central events than for peripheral. 

In the region (ii) the exponent in (j!9|) can be neglected which yields 

^- = S A xG{x u Q 2 ) ' lee . (24) 
dy (2vr) 2 a em 

With the help of ()22|) we find the behavior of the nuclear modification factor at high ener- 
gies/forward rapidities 

M * <<Q " (25) 

where we the subscripts A and p are introduced to distinguish Q in nucleus and proton. It 
gets suppressed both as a function of energy/rapidity and centrality. 

C. Limiting fragmentation of J/ty and hidden parton scaling 

1. Total cross section 

As has been mentioned in the Introduction a naive collinear factorization approach implies 
that the total J/^> production cross section (as well as the cross section of any other hard 
process) is proportional to the product of parton distribution functions of the proton xf(x{) 
and of the nucleus a; /a (#2)- Moreover, if the coherent effects are neglected, then x/a^) = 
Axf(x2), i- e. the nuclear effect factorizes out. Therefore, in the collinear factorization the 
total cross section is proportional to 

da(y) 



dy 



ex xf(Q<?/y/Z)xf{Qe-v/y/Z), (26) 

coll 
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where Q is a typical scale of the hard process. At some other energy yfs* the cross section 
is: 

oc xf(Qey + Y /V^)xf(Qe~y + Y /V^), (27) 

ay coil 

where Y = (1/2) ln(s'/s). However, such dependence on energy contradicts experimental 
data on inclusive particle production in which do~(y)/dy exhibits scaling with Y + y in 
forward region y > 0. Analogous phenomenon for soft processes is known as the "limiting 
fragmentation" . 

The scaling of da(y)/dy with Y + y or, equivalently, with x\ is a natural consequence 
of saturation of the nuclear wave function at y > 0. It means that the cross section has 
the same shape as a function of x\ at different energies, see Fig. Efb). In the case of J 
production it is manifest in Eq. (|24j). see Fig. Efb). We can see that wee partons of the 
nuclear wave function are saturated and hence do not contribute to the fragmentation in 
the forward rapidity region. The parton scaling in the nucleus is effectively hidden at small 
X2\ we thus use the term hidden parton scaling to describe the universal scaling in X\. Since 
xp = x\ when x<i <C 1 the hidden parton scaling of da(y)/dy is equivalent to xf scaling in 
the same kinematical region. 



2. Nuclear modification factor 

As we noted above, the A-dependence of the cross section in the collinear factorization 
approach is trivial. Consequently, the nuclear modification factor (}2*Tj) equals unity; it thus 
scales with x 2 and x\. Coherent scattering of the proton in the nucleus (manifesting itself in 
the dependence of the scale Q on A) breaks this scaling. This can be seen in Eq. ([25)1 . which 
shows an explicit dependence on x\ and x%- Note that the difference between xG(xi,Q\) 
and xG(x\, Q^) is largest in the proton fragmentation region x x — ► 1 corresponding to very 
low xo. Here we expect the strongest violation of x<i scaling in agreement with experimental 
data [56(. 

If we compare R p a{J/^) at close energies we find an approximate xf scaling. This scaling 
originates in the slow dependence of the gluon distribution on energy, Q 2 ~ xG ~ s A / 2 ; as 
a result, both x\ and xp scalings approximately hold for close energies, see Fig. El (The 
scaling is broken at very forward rapidities x\ — > 1, where the sensitivity to the variation of 
the scale Q is enhanced by the power fall-off of the parton distributions, (1 — xi) n .) We thus 
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explain the Xp scaling observed for the SPS and Fermilab fixed target energies. Variation 
of xG between SPS and RHIC energies produces a much stronger violation of xp scaling as 
seen in Fig. |H1 

In our discussion so far we have always assumed that the proton wave function never 
saturates. However, the gluon saturation in the proton has been likely observed at HERA 
(see e. g. Ref. at x's somewhat smaller than those accessible at RHIC. However, at the 
LHC much smaller values of 22 can be reached at both central and forward rapidities and 
the proton wave function may form the Color Glass Condensate. In that case R p a becomes 
a function of only x\. This statement is general for all hard process, the J/ty production 
being a particular example. We thus predict that the hidden parton scaling, or the scaling 
in x±, will become a universal feature at the LHC energies. 



V. COMPARISON WITH RHIC DATA 



A. A model 



We will now compare our calculations with the experimental data using Eq. (fTT)|) which 
is the total J/^> production cross section in the quasi-classical approximation. To take into 
account the quantum evolution effects on the scattering amplitude we use a model suggested 
in [5^]. It gives a good description of inclusive particle production at RHIC. In that model 
the quantum effects in the scattering amplitude are parametrized as 

1 - e" n2 - 1 - e -^ Q2) (28) 

where f2 = \ Q?,(x>?) r 2 , and ■y(y,Q 2 ) is the anomalous dimension. Its explicit expression 
can be found in Ref. [5J|. Let us only note here that 7 is chosen in such a way as to satisfy 
the large Q, fixed y as well as large y, fixed Q asymptotic of the DGLAP and the BFKL 
equations up to the NNLO terms. 

Let us now list important assumptions which we made in deriving ()19j). First, the scat- 
tering amplitude is calculated in the large N c approximation. Second, in our calculations we 
assumed that the J/*f> wave function is non-relativistic. In this approximation the g* — J '/\& 
overlap function takes a simple form shown in (|14)1. Relativistic corrections due to Fermi 
motion strongly depend on the charmed quark mass. On other hand in the case of diffractive 
j/\[r production in DIS it was observed that these corrections are almost energy independent 
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Hence, we consider the effect of Fermi motion as an uncertainty of the wave function 
normalization 5^ K F . 

Third, in derivation of (fE^) we assumed that only c and c interact with the nucleus. 
Of course the valence quark and the gluon can interact as well. These processes give a 
contribution of the same order as the one we discussed in the previous subsection, and their 
parametric dependence on energy and atomic number is the same. Therefore, we can also 
take them into account in the overall normalization factor K F (3]. 

Finally, there are parametrically small corrections to Eq. (jl9j) due to contributions of 
the real part of the amplitude and off-diagonal matrix elements. These corrections are 
numerically insignificant and we have neglected them. 



B. Attenuation in a cold nuclear matter 

We have argued in Sec. Ill Bl that our result (fTT?jl is applicable at rapidities y > at 
RHIC. If we would like to compare it with the lower energy data we need to include the 
effect of absorption in the cold nuclear matter due to finite coherence length. In that case 
the produced cc pair and later J/\l/ itself can inelastically interact with the nuclear matter. 
This reduces the cross section (fT9j) by a factor (see e.g. [30]) 

= e -^^W, (29) 

where a^ N is the inelastic cc ( J/\l/)-nucleon cross section, p = 0.17 fm~ 3 is the nuclear den- 
sity and L(b) is the length traversed by J/\l/ in nuclear matter at a given impact parameter 
b. The nuclear absorption factor S$ has been studied in detail in the framework of Glauber 
approach Ref. We normalize in (J2~9~j) accordingly. 



C. Results of numerical calculations 

We have performed a numerical calculation using the model described in Sec. | V A | The 



47, m. The 



parameter k has been fixed at k = 0.2 GeV 2 , consistent with the analysis of 
overall normalization factor K F in the production cross section has been fitted to the RHIC 
data from 0, IsfJ. 

The result of numerical calculations of the total J/\l/ cross section as a function of rapidity 
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FIG. 4: Total inelastic cross section for J/^f production (|19|) as a function of (a) rap idity y, (b) 
shifted rapidity y + Y, where Y = (1/2) ln(s'/s). Data points are taken from Ref. 561 ] . 

is shown in Fig. HJa). We observe a reasonable agreement with the PHENIX preliminary 
data 

In Fig. IHb) we show the "hidden parton scaling" due to the saturation of the nuclear 
wave function, see Sec. II V CI We observe an effect similar to the "limiting fragmentation" 
in the total inclusive cross section. Let us emphasize that such a "hidden parton scaling" is 
a general feature of all hard processes at high energy in the saturation picture. It is due to 
the saturation of partons in the nucleus. 

In Fig. we show the nuclear modification factor as a function of rapidity and centrality. 
In agreement with our discussion in Sec. II VI B t a( J/^/) is suppressed at forward rapidities 
and enhanced at midrapidity. The functional dependence on rapidity and centrality isgiven 
by Eqs. f|23 |h (|25 |) . Our results are in a qualitative agreement with the PHENIX data 

HQ- 

Finally, in Fig. |B] we present our calculation of the exponent a defined as follows: o v a = 
A a a pp . As discussed in Sec. II V CI we expect an approximate xp scaling at close energies due 
to the slow dependence of the saturation scale on energy. At smaller xf nuclear absorption 
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FIG. 5: Nuclear modification factor Rj/^ as a function of (a) rapidity y, (b) iV^jj at different 
rapidities. Solid lines: numerical calculation. Data points are from Ref. 56]. 

plays a significant role; because of formation time effects, it suppresses the J/^f production 
at SPS stronger than at Fermilab and contributes towards improving the xf scaling. In 
the same figure we plot our prediction for LHC At such high energies the proton's wave 
function saturates, which results in the exact X\ scaling, as discussed in Sec. II V CI Thus, in 
the region x\ ^> x 2 we predict Xp scaling of a at energies >5 TeV. 

We would like to draw the reader's attention to the fact that a ~ 2/3 at the highest 
energy in Fig. This is just the value of a one expects to measure in the collision of two 
black disks: one of radius R p another of radius Ra- At LHC energies the absorption starts 
only at y < —3. Therefore, from our discussion in Sec. HVBl we expect that a will approach 
4/3 at y < 0. This behavior will be best seen in a plot of a as a function of y since the plot 
versus Xf exponentially expands the fragmentation region, while shrinking the interesting 
central rapidity one. 



VI. SUMMARY AND CONCLUSIONS 



In this paper we have analyzed the J/^> production in p(d)A collisions at high ener- 
gies. We have pointed out that due to the large coherence length associated with the J/^f 
production and formation the dominant production mechanism in the case of heavy nuclei 
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FIG. 6: a as a function of xp for different energies, a was extracted by fitting function A a to the 
calculated o v aIo vv . Bands show the uncertainty of the fit. 

is a double-gluon exchange which leads to a significant enhancement of the cross section 
at backward rapidities. In the case of DIS this effect has been discussed in Ref. |35fl. and 
the relevant arguments for hard processes at forward rapidities have been given in [8j. At 
forward rapidities the J/^> production is suppressed in much the same way as open charm 
is (scj. This is due to saturation of gluons in the nuclear wave function, see Fig. |3Jb) . The 
transition between the enhancement and the suppression regimes in J/\l/ production hap- 
pens when Q s « M^. Eq. (|T9~j) provides an analytical formula for the total cross section of 
J/\[r production in the region where Tp > Ra- 

The dominance of the double-gluon exchange mechanism for the case of heavy nuclei 
has an important implication. It means that J/^f is produced directly in the color singlet 
state. A clean way to measure the contribution of the double gluon exchange from different 
nucleons is to measure J/ty production in proton-deuteron collisions and to trigger on 
the final fragmentation state of the deuteron. A signature of the double-gluon exchange 
mechanism will be the absence of intact nucleons j57]. Such an experimental study can shed 
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further light on the problem of J/^f production. 

We have argued that the total J/*f> production cross section exhibits the "hidden parton 
scaling", see Fig. Efb). Owing to the saturation of gluons in the nucleus, the cross section 
becomes independent of x<i- Thus it scales with x±, or, equivalently, with y+Y = lnxi+const. 
In other words, da{xi) / dx\ has a universal shape for different energies. 

In Sec. II V CI we studied the xp scaling phenomenon observed at lower energies. We found 
that this scaling holds only for close energies since the scaling violating factor is a slow 
function of energy: ~ s 25 . However, at energies as high as LHC energy we expect the 
saturation of gluons in the proton which will manifests itself in the effective disappearance 
of the factor s 0,25 and will result in the exact x\ scaling. In turn, an exact x\ scaling can be 
considered clS cL SI gnature of the saturation in the proton. 
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